Effect of Random On-Site Energies on the Critical Temperature of a Lattice Bose Gas 
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, We study the effect of random on-site energies on the critical temperature of a non-interacting 

Bose gas on a lattice. In our derivation the on-site energies are distributed according a Gaussian 
' probability distribution function having vanishing average and variance v^. By using the replicated 

action obtained by averaging on the disorder, we perform a perturbative expansion for the Green 
functions of the disordered system. We evaluate the shift of the chemical potential induced by the 
\^ ' disorder and we compute, for v1 <^ 1, the critical temperature for condensation. We find that, for 

^ large filling, disorder slightly enhances the critical temperature for condensation. 
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I. INTRODUCTION 



The study of the effect of disordered potentials on quantum particles is, since decades, an important field of 
research, one of the main motivations being the understanding of how electrons in disordered systems localize P, 
. The investigation of disorder effects on bosonic systems, motivated by experiments of adsorption of ''He in 

' porous media has been also very active, especially in relation to the issue of understanding how superfluidity 

^ ; ' properties are modified in random environments @, @, 0]- In this respect, ultracold bosonic trapped gases Q 
' provide a good experimental setup to study the effects of disorder on bosonic systems. In an ultracold Bose gas, 
^ \ disorder can be induced both by a laser speckle [13, Eli or by an incommensurate bichromatic potential [13] (i.e., an 
I auxiliary, incommensurate, lattice added to an optical lattice). These techniques recently allowed for the experimental 
■ O observation of Anderson localization for matter waves in a random potential [l^, [l^ . 

An important resource to control the properties of ultracold bosons is provided by the possibility of superimposing 
on them optical lattices [T5| . allowing for a fine tuning of the ratio between kinetic and interaction energies. The 
effective Hamiltonian for bosons in a deep optical lattice is the Bose-Hubbard Hamiltonian [l,[l3], i-e., a tight-binding 
model characterized by a kinetic term describing the hopping of bosons between neighbouring sites of the lattice 
and an interaction term proportional to the s-wave scattering length between bosons. In absence of disorder, the 
Bose-Hubbard Hamiltonian has a superfluid (Mott insulator) ground-state for large (small) ratio between kinetic 

^ ^ and interaction energies. Random on-site disorder together with a strong interparticle interaction between lattice 

f — . bosons is expected to induce a Bo se g lass phase [^ [l7| which has been recently studied both from a theoretical 



Tj- ■ [S [il, M lli 122, l23|, IM 125|, l26|, [27|, HE HI and experimental [HI side. In the opposite limit of very weak interaction 
\ at T = the Anderson localization has been studied [3] '■ the study of this limit is feasible thanks to the possibility of 
' tuning the s-wave scattering length a of a ^^Rb Bose gas with high precision, and setting it almost to zero [s^l- When 

^\ ' an optical lattice is superimposed to an ideal Bose gas, for large values of the laser power, the system is described by 
■ a boson-hopping model, i.e., the Bose-Hubbard model [EHBl without the interaction term. 

' In [3l| the thermodynamic properties of a Bose gas in a disordered lattice have been analyzed and the shift of the 
critical temperature for condensation induced by the disorder was determined. For a three-dimensional lattice it has 
been found that the shift of the critical temperature depends on the filling, i.e., on the number of particles per lattice 



^ , site. When disorder affects the hopping rates between neighbouring sites, the critical temperature Tc is enhanced for 
- - ■ large filling / and the shift does not sensibly depend on / |3l|; at variance, for small /, decreases. In presence of 
random on-site energies, increases for large filling, but much less than for for bond-disordered lattices, resulting in 
very small shift of Tc for small disorder. These results should be compared with the findings for a continuous (i.e.. 



34 , |35|, |36| : without any confining 
H. In Refs. [li,!!^ the properties, 



without optical lattice) Bose gas in presence of a disordering potential [3j, |3S 
potential, it has been shown that the critical temperature decreases with disorder 
both at T = and at finite temperature, of a continuous Bose gas in a random potential were investigated in the 
frame of a self-consistent mean-field approach. The shift 5Tc of the critical temperature has been computed in a 
Hartree-Fock approach by Lopatin and Vinokur [sl] as 



2 



where p is the density, m the mass of the bosonic atoms and Tc the critical temperature for condensat ion i n absence 
of disorder. In Eq. ([1]) /C is the strength of the disorder: the disorder potential u{r) has averages u{r) — and 
■it(r, F) = K.5{r — f) (the bar denotes the average all disorder configurations). The result ([T]) has been confirmed in 
[33| by one-loop Wilson renormalization-group calculations. The shift of the critical temperature for the continuous 
Bose gas has been also recently computed in [sl] using the Popov method, obtaining a value for the relative shift 
STJT^°^ which differs by a factor of 1/2 from the resuh A discussion of the reasons for such difference is contained 
in [33 • In presence of an optical lattice, the decrease of the critical temperature is found for small filling, for which 
the result ([T]) is retrieved [3l|: however, on the lattice, different behaviours are possible as a result of the interplay 
between discreteness and disorder [3l|. 

The aim of this paper is to provide details of the results presented in [31], focusing on the physically relevant case 
of random on-site energies. We shall present a detailed derivation of the computation of the Green functions at the 
first order in the disorder parameter, and discuss the contributions from higher order terms. The plan of the paper is 
the following: in Section II the model Hamiltonian is presented and the results in absence of disorder are reviewed. 
In Section III the computation of the replicated action is presented. In Section IV, the shift of the chemical potential 
induced by random on-site energies is computed by a perturbative expansion of the Green functions and the resulting 
shift of the critical temperature is determined and discussed. Our conclusions are in Section V. 



II. THE MODEL HAMILTONIAN 

When there are random on-site energies, non-interacting bosons on a three-dimensional lattice are described by the 
Hamiltonian 

^ = E {^rbj + blk) + E ^^b.- (2) 

In Eq. i is the tunneling rate between neighboring sites, the lattice sites are denoted by the operator hi (hi) 
destroys (creates) a boson in the lattice site i and the sum is on all the distinct pairs of neighbouring sites. The 
number of sites is denoted by Ns'- for a cubic lattice with of linear size L, Ns — ■ The total number of particles is 
Nt and the filling (i.e. the average number of particles per site) is given by 

In the Hamiltonian ([2]), the random on-site energies are accounted for by the e^'s which are regarded as random 
variables with vanishing average and variance v^t^ . We shall consider a small disorder {v^ ^1) and we assume that, 
at each site, the disorder has probability distribution P{ti) given by 

P[e.)^^^e-^y^<. (4) 

When there is no disorder (e^ = 0), the Hamiltonian 12]) reads 

Hi^)^~tY,{b%+b\K) (5) 

(hereafter, the superscript '^"^ is used to describe bosons when on-site disorder is absent). The critical temperature 
for condensation Ti'^^ can be determined in the usual way Q : one introduces the chemical potential p, to enforce 
the conservation of the total number of particles and replaces //f"' with K^'^^ = ij'"' — /iiV, where N — SJS; is 

the total number operator. Requiring {N) = Nt and performing the thermodynamical limit {Ns,Nt — > 00 at fixed 
filhng /), one gets 

1 ^ 

BZ (27r)3 e0iE^-^^) - 1 ^ ' 

where the integral is on the first Brillouin zone and 

i^k = —2t (cos fca: + COS ky + COS k^) (7) 
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are the single particle energies (periodic boundary conditions have been assumed). The critical temperature ri"-* is 
then determined by requiring //(Ti"-*) = /^^ = — 6t, which yields 



(where f}^"' - l/fcsTrO 
In [m it is shown tha 
lowest order of the Taylor expansion of the exponential in Eq. ([8]). One gets 



In [31[ it is shown that a reasonable estimate of Ti"', valid for large filling, may be obtained by keeping only the 



-W(l)' ^' 
where W(l) ~ 1.516386 is the value of the three-dimensional generalized Watson's integral 

W(z)^ / -4 (10) 

Jbz (2^)' 1-^ELiCosfc, 

computed in z = 1 [11]. In [3l| the value of Ti^"* given by Eq. ([9]) is compared with the numerical solution of Eq. ([8]), 
showing that this estimate is quite good also for intermediate values of the filling: e.g., for / > 5, the relative error is 
less than 5%. 

III. EFFECTIVE REPLICATED ACTION 

For the Hamiltonian ^ with random on-site energies, the partition function can be written as a path integral (39j 

Z = f l[V^,{T)V^*{T)e-'' (11) 

(h — ks ~ 1), where the action S is given by 

S^J dri^Y^ip* (^-^^ fi^ ip, + Y^eiip*ipi-tY^ (12) 

Of course, the partition function (jlip depends on the on-site energies distribution, i.e., Z ~ Z{{e}). The fields ip's 
describe Bose particles. 

As it is usual when one deals with a disordered system, one introduces N replicas of the system and, after averaging 
the action on the disorder, performs the limit iV ^ [40]. This allows to write down a perturbative expansion of the 
Green functions. Labeling the N replicas by a = 1, • • • ,N, one has 



■ (13) 

The averaged effective partition function is given by 

Z^WV) = / n de,P{e,)Z^{{e}l (14) 
where P{ei) is the Gaussian probability distribution After integrating over the on-site energies e^'s, one gets 

'W{{c\) = I llV^nr)Viprir)e-'^^f, (15) 
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where the effective rephcated action Seff is given by 




As one can see from ()16p . in the rephcated action the disorder enters as an effective attractive interaction between 
rephcas with strength given by vlt^ /2. 



IV. SHIFT OF THE CRITICAL TEMPERATURE 

In this Section we determine the shift 5Tc = Tc ~ Tc^^ of the critical temperature for non-interacting bosons in a 
three-dimensional lattice with random on-site energies. To do this, we use the replicated action (jl6p to write down a 
perturbative expansion in v'^ of the Green functions of the disordered system: this allows to compute both the shift 
of the chemical potential induced by the disorder and the shift 6Tc. 

The relevant Green functions are 

Gij{z) ^ {(pi(p*)v^=io (18) 

where z — fi + iojn and ujn are the bosonic Matsubara frequencies. With Q (G) we denote the Green function when 
there is not (there is) disorder. Upon introducing the self-energy function 51 (k, z), one may write the Dyson equation 
in Fourier space as 

G-i (k, z) = (k, z) + [5:(k, z) - 5iA , (19) 
where the Fourier transform of Q is given by 

g-i(k,z) = £;k-z (20) 
and (the shift of the chemical potential) is defined as 

= 5:(k = 0,0 = ^). (21) 

Bose-Einstein condensation occurs when 

G(k,z)-i = -6t- Ai+ [5:(0,m) - M = (22) 

k— 0,2— 

which, at zeroth order in v^, is solved by 

= -6t : 

thus, at any order, the chemical potential at is given by /ic + 5^^ = —Gt + Z(— 6t). 

To compute Eq. (|18p at the first order in u^, one needs to evaluate the first terms of the Taylor expansion of e~^''ff . 
Using Wick's theorem one gets 

The quantities 

are diagonal in the replica indices. From Eq. (j23p one gets 

Glrn = Ghn + ^(1 + N) Qu g^^ G „n + 0{vt), (24) 
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which for ^ leads to 



Glm = Qlm + Qu Qu Gim + 0{v'^). 



(25) 



From Eq. (|25p and using Eq. ^Tl} one gets for the self-energy 
Ii(k, z) = Zi(z) = -vlt^Gu = -vy 



dk' 1 



where the subscript i denotes that we are including contributions up to the order u^; furthermore 



(26) 



(27) 



with \N{z) defined in Eq. pU|) . From Eq. one may readily compute at the critical point the shift of the chemical 
potential 5^ci at the first order in v^: one finds 



0.25 f ft. 



(28) 



As one can see from Eq. (1^^ . due to the contribution given by the self-energy to the chemical potential, the filling 
fraction / is changed in presence of disorder to / -f 5fi, with 5fi given by 



5h = lim T^Ve— / ( ^— — 

^ Jbz (27r 



Z — llc-\-'>'^n 



r)3 {E,,-zf 



Substituting 



BZ 



(2^ (-Bk - ~ 



in Eq. ([29)l . one gets 



= hm Tc 



F(^)-F(/x,) -F(z) 

J oz 



Z — flc.-\-iuJn 



(29) 



(30) 



(31) 



For large values of the filling /, the dominant contribution in Eq. (pij) is given by the lowest Matsubara frequency, 
which yields 



S.h = vi Te lim <; F(z) - F(/i,)l |-F(z) 

z^Mc I L 1 oz 



(32) 



For an explicit computation of (|32p . it is most useful to define the dimensionless variable z = z/t and the function 
F(i) = -(l/5)W(-S/6). Since ii^/t = -6, one has 



where the numerical coefficient A is given by 



A = lim 



F(i)-F(-6) |.F(5) 



(33) 



(34) 



In order to compute A, one may use the expressions of F and W in terms of elliptic integrals given in [4l| : it is 

2 ' 



W(i) = - 



(35) 
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Fi~z) = W(l) -^Jl-^\+ 0{r,^) (37) 



and 



Inserting Eqs. jST])-!!!!]) in Eq. ((341), one readily finds 



As a result of Eq. ([9]), one notices that, for large /, there is a linear relation between Tc^^ and /, which is given by 



^(^)t(0)./. (40) 



6t 



At the order v^, taking into account the shift in the particle density given by Eq. (|39p . one has 

W(l) 

m 

From Eqs. and (HT1) and using Eq. ([M)) . one obtains 



Tc + Sfi = /. (41) 



Tc = T!:°Ul + -r4n77Tvl)^ (42) 



from which 



167r2W(l) V ' 



- 0.0125 (43) 



167r2W(l) 



To go beyond the first order results, one may be tempted to use the approach based on the momentum space 
renormalization of the vertex function: within ladder approximation, the renormalized vertex V(k, zi,Z2) can be 



written for zi = Z2 = z as 



where the polarization function is given by 

n(k,z)==/ -j^- ^. (45) 



As a result, the self-energy II(z) is given by 
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and it contributes to the equation for the filhng fraction as 



/= hm_r,^e-*""^ 



BZ 



(2^ 



Z(z) - (5^c 



(47) 



At the first order in u^, the vertex function (|44p is simply given by v^t^ and the corresponding self-energy is '^{z) — 
—v1t^F{z), as given by ((26l) : using Eq. (|47)) . one, of course, retrieves for the shift of the critical temperature the result 
()43|) . If one attempts to compute the contribution to / coming from higher powers in vq, one is immediately faced 
with problems arising from infrared divergences. In fact. 



n(k = 0,z) 



9F(z) 
dz 



diverges at z = ^c- The situation is similar to what happens for the computation of the shift of the critical temperature 
due to weak rep ulsive interactions, where the perturbation theory also fails to give a finite expression for (STc, as 
discussed in I43l. 



V. CONCLUSIONS 



In this paper we investigated the effect of random on-site energies on the critical temperature of a non-interacting 
Bose gas on a lattice. By using the replicated action obtained by averaging on the disorder, we performed a per- 
turbative expansion for the Green functions of the disordered system. We computed, for <C 1, the shift of the 
chemical potential induced by the disorder and the critical temperature for condensation. For large filling, the critical 
temperature is slightly enhanced with respect to the situation in absence of disorder: the relative shift STc/Tc^^ is 
given by STc = aUg, where the numerical coefficient a is positive and small (a w 0.0125). 
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